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Abstract 

The  plasma  shape  and  location  of  a  sharp  boundary 
belt  pinch  is  determined  by  the  dimensions  of  the  sur- 
rounding conducting  shell  and  by  thermodynamic  param- 
eters.  It  is  shown  that,  except  for  simple  scaling, 
there  exists  a  one— parameter  family  of  equilibria.   The 
plasma  shape  is  uniquely  determined  by  initial  parameters 
and  wall  current.   The  variation  of  the  plasma  shape 
during  an  adiabatic  compression  is  calculated  as  are 
plasma  current,  temperature,  field  strength,  etc.   As  the 
compression  ratio  of  the  plasma  is  Increased  beyond  two, 
the  ratio  of  plasma  height/width  decreases  to  unity. 
For  high  compression  ratios,  the  column  cross— section 
becomes  circular.   Explicit  solutions  are  obtained  by 
conformal  mapping. 
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I .   Introduction 

The  desire  for  stable  toroidal  plasma  equilibria  has  led 
to  studies  of  plasma  configurations  with  non-circular  toroidal 
cross-sections.    Many  of  the  efforts  to  determine  equilibria 
for  such  configurations  have  been  numerical  investigations 
which  assume  a  sharp  boundary  toroidal  plasma,  supported  in 
equilibrium  by  a  vacuum  magnetic  field.   These  calculations, 
which  consider  a  variety  of  plasma  current  density  distributions 
offer  some  surprising  results.   For  example,  plasma  compression 

ratio  limitations  of  the  order  of  two,  for  highly  non-circular 

2  ^ 
plasmas,  have  been  observed.  '    This  limitation  may  be 

attributed  to  the  lack  of  circular  symmetry  of  the  plasma  cross- 
section.   In  determining  such  limitations,  the  aspect  ratio  is 
less  important  than  the  ratio  of  dimensions  of  the  cross- 
section  of  the  outer  boundary  of  the  torus. 

In  this  paper,  we  present  an  analysis  of  configurations 
in  which  all  the  current  in  a  uniform  plasma  column  flows  in 
the  interface  between  the  plasma  column  and  the  vacuum  magnetic 
field  which  supports  it  in  equilibrium.   Numerical  studies 
Indicate  that  this  configuration  offers  maximal  compression 
ratios  for  a  belt  pinch  column.   We  simplify  our  problem 
by  neglecting  toroidal  curvature  effects.   That  is,  we 
replace  the  toroidal  configuration,  assumed  to  have  a  perfectly 
conducting  closed  toroidal  shell  of  rectangular  cross-section 
(see  Flg.l),  by  a  straight  configuration.   The  toroidal  axis 


of  symmetry  is  replaced  by  a  plane  of  symmetry.   The  closed 
toroidal  outer  boundary  is  replaced  by  two  perfectly  conducting 
straight  sections  symmetrically  located  on  opposite  sides  of 
the  plane  of  symmetry  (see  Fig.  2).   These  straight  sections 
carry  oppositely  directed  "toroidal"  current. 

Toroidal  current  in  the  outer  boundary  induces  a  toroidal 
sheath  current  at  the  plasma  boundary  v/hich  in  turn  induces 
a  poloidal  vacuum  field  in  the  region  between  the  plasma 
column  and  the  outer  boundary.   An  increase  in  this  field  alters 
the  pressure  and  shape  of  the  column  in  order  to  provide  local 
pressure  balance.   At  the  interface 


2Tr  ^  col  -  2Tr  ^  vac         ■       ^^^ 


is  t h e  strength  of  the  vacuummagnetic  field  at  the  plasma- 
vacuum  interface.   The  fundamental  problem  is  to  determine  the 
shape  and  location  of  the  interface.   An  explicit  determination 
is  obtained,  in  terms  of  elliptic  integrals,  using  conformal 
mappings.   This  standard  technique,  due  to  Kirchoff   for  the 
solution  of  free  boundary  problems,  was  first  applied  to  plasma 
containment  in  the  study  of  cusped  geometries,   and,  more 
recently,  in  the  determination  of  the  location  and  shape  of  the 
plasma  column  during  an  adiabatic  tokomak  compression.    For 
boundaries  which  are  very  high  compared  to  their  width  (o  >>  1, 
see  Fig.  2),  we  show,  using  an  asymptotic  analysis,  that 


If  the  ratio,  column  height/column  width  Is  to  be  large, 
then  the  toroidal  current  used  to  compress  the  column  can  not 
be  too  great.   Very  high  compression  can  only  be  achieved  by 
a  sacrifice  of  column  elongation. 

Numerically,  we  find  that  an  equilibrium  plasma  column 
with  a  compression  ratio  less  than  two  has  dimensions  comparable 
to  those  of  the  confining  outer  boundary.   As  the  compression 
ratio  is  increased  beyond  two,  the  ratio  of  plasma  height/ 
width  decreases  to  unity.   For  high  compression  ratios,  the 
column  achieves  an  essentially  circular  cross-section. 
Detailed  results  are  presented  In  Section  IV  for  ^  =  5-9,  the 
dimensions  of  the  belt  pinch  at  the  Max  Plank  Institut  fur 
Plasmaphysik  at  Garching  bei  Munchen. 


II .   The  Free  Boundary  Problem 

In  the  vacuum  field  region,  the  poloidal  field  satisfies 
the  equations 

V  X  B  =  0  (2) 

V  •  B  =  0     .  (3) 

These  equations  imply  the  existence  of  a  scalar,  $,  not 
necessarily  single-valued,  such  that  in  the  vacuum  region 

B  =  -V$  (4) 


(5) 


In  view  of  the  symmetry,  it  suffices  to  determine  $  in  the 
vacuum  region  contained  in  S  =  {(x,y)|0  <x<H,  0<y<R} 
(see  Pig.  3). 


For  convenience,  we  introduce  normalized  variables.   We 
normalize  the  poloidal  field  so  that  its  strength  is  unity  at 
the  plasma  boundary.   We  define  iJj   to  be  the 
negative  of  the  ratio  of  the  vacuum  poloidal  flux  to  the 
poloidal  field  at  the  plasma-vacuum  interface.   The  normalized 
poloidal  field  will  be  denoted  by  b.   Its  x  and  y  components 
will  be  denoted  u  and  v  respectively.   The  normalized  potential 
is  (j)(x,y)  and  its  harmonic  conjugate  (the  stream  function) 
is  tj;(x,y  )  .   Thus 

b  =  -V())  (6) 
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V  =  -t),^  =  +^^  (8) 

A(t)  =  Aijj  =  0  .  (9) 

The  boundary  conditions  which  (}>  must  satisfy  are:   (i)  on  the 
straight  line  segments  through  y  =  R  and  x  =  H,  ^  =  0  and 
■r-^  =  0  respectively;   (ii)  on  each  of  the  remaining  straight 
line  segments  bounding  the  vacuum  field  region  contained  in 
S,  (p    is  constant;   (ill)  at  the  interface  |  V(|)  |  =  1. 

It  happens  that  for  our  normalized  potential,  there  is  a 


method  used  to  determine  (J),  and,  at  the  same  time,  the  free 
boundary,  is  due  to  Kirchoff  and  utilizes  conformal  mapping. 
Consequently,  we  introduce  the  complex  variables  z,  w,  and  f 
defined  by 

z  =  X  +  iy  (10) 

w  =  u  -  iv  (11) 

f  =  (})  +  14)      .  (12) 

Assuming  the  real  and  imaginary  parts  of  w  and  f  to  be  twice 
continuously  dif f erentiable  functions  of  x  and  y,  we  find  w  and 
f  to  be  analytic  functions  of  z.  Moreover,  from  Eqs.  (10)-(12) 

w(z)=-|f        .  (13) 


parametrically  by 


z  =  -  1   -  df  +  1A-,  (lU) 

J   w        1 

0 

f'*°   1 

iR  =  -       -  df  +  lA^    .  (15) 

I      w        1 


0 
If  z  Is  on  the  free  boundary,  we  have 


(16) 


0 


Thus,  the  free  boundary  is  determined  If  w  can  be  determined  as 
an  analytic  function  of  f.   This  is  accomplished  by  a  conformal 
transformation  which  maps  the  range  of  f  onto  the  range  of  w. 
These  two  ranges  are  shown  in  Figs.  (4)  and  (5).   The  mapping 
of  the  range  of  f  onto  that  of  w  must  preserve  the  correspondence 
of  points  of  common  origin.   Thus,  using  Figs.  (3),  (^)    and  (5), 
we  obtain  the  following  table  for  the  mapping.   In  this  table 
we  have  defined  the  number  [(})]  as 


(17) 


where  the  Integral  is  taken  so  that  the  vacuum  field  region 


satisfy  0  <  a.  <  1. 


z 

w 

f 

iR 

-"l 

i^ 

lA^ 

-1 

0 

h 

-i 

J[*] 

H 

-ia2 

^Ul+i'l'o 

H+iR 

0 

-a^  *  ^*o 

TABLE  I 


This  table  determines  a  unique  analytic  transformation  taking 
the  Indicated  domain  of  f  onto  that  of  w.   In  particular. 
Introducing  the  auxiliary  complex  variable  X      It  Is  possible 
to  write  w  and  f  parametrlcally  as  the  Image  of  the  half  plane 
Im  A  >  0 


(18) 


(19) 


In  equation  (19)  the  quantities  A   and  X   are  the  values  of  A 
associated  with  z  =  IR  and  z  =  H  respectively  (I.e.  f  =  lip 
and  f  =  ■n-[<t>]  respectively).   The  numbers  A  =  0  and  A  =  1  have 
been  assigned  to  z  =  lA-j^  and  z  =  Ag  respectively  (see  Fig.  6). 
The  corner   z  =  H+IR  has  been  mapped  to  infinity  in  the  A  plane 


The  quantity  r  is  determined  by  the  condition 

1^^  =  I   If  dX        .  (20) 

0 

In  equations  (l8)  and  (19),  square  roots  are  taken  with  branch 
cuts  along  the  negative  real  axis  and  principle  value  along  the 
positive  real  axis. 

Equation  (18),  which  maps  the  range  of  w  onto  the  half- 
plane,  may  be  obtained  by  first  mapping  the  w-domain  onto  the 
first  quadrant  with  a  rotation,  followed  by  a  squaring  and  a 
fractional  linear  transformation.   A  second  squaring  covers  the 
half-plane,  and  a  second  fractional  linear  transformation 
assigns  points  as  described  above.   Equation  (19),  mapping  f 
onto  the  upper  half-plane,  is  obtained  using  a  Schwarz- 
Chrlstoffel  transformation. 

Equations  (15),  (l6),  (l8)  and  (19)  give  the  shape  of  the 
free  boundary  explicitly  as  an  elliptic  integral.   Integrals  of 
this  type  are  readily  evaluated  using  standard  numerical 
proceedures .   At  critical  parameters  ranges,  where  numerical 
calculations  are  awkward,  asymptotic  approximations  may  be 
found. 


III.  Asymptotics 

After  changes  of  variables  in  equations  (l4),  (l8)  and 

(19),  we  obtain,  writing  e  =  -A„  and  1 

n 


0) 


/2- 


d? 


_1  AI?  72+^(1+?)  I  72(1+6)    y2(l  +  6)-(l+OJ 

dg      1  1 


.^  A-5^  ^1+|(1+C)  ^2{l+&)  +  c{l  +  0 


(21) 


/2- 


d5 


(e/1+6) 


1/2 


_^  ^l-C^  -^1  +  5+ff^  I  >^2(l  +  6)    v/   l+.&+2e 
+  1    dl       1 


L-r^  A+^(l  +  E)  ^S 


1-C    ^1+^(1+^)  ^2(l+6)  +  e(l  +  e) 


(22) 


d? 


-1  vi-g  yi+g+2£  V2(i+5)-(i+iy 

•  +  1    dE,  1  1 


h-e  A+^ 


(23) 


+  1 


dC 


-1/T+I  •l+g+2£  /2(l+6)-(l+gT 
1    dg        1 1 


_^  A-g^  A+|(l+g)  ^2(l  +  6)  +  e(l+e) 


(24) 


We  shall  assume  that  the  dimensions  of  the  outer  boundary,  H 
and  R,  have  been  given.   Then  equations  (21)  and  (22)  determine 
permissible  values  of  e  and  6  (i.e.  X   and  A„)  as  functions  of 
■<\)    .      Thus,  for  ^      given,  equations  (l4),  (18)  and  ( 19 )'  determine 


the  free  boundary.   In  particular,  equations  (23)  and  (24)  give 

the  locations  of  the  extremities  of  the  free  boundary.   We 

H 


In  the  limit 


ir—T  <<    1.   The  condition  t—-t-   <<  1  is  a  necessary 

1+0  1+0 


equations  (21)  and  (22),  we  find 


+  1 


d? 


•^^TT 


1  +  6 


(e/1+6) 


1/2 


A-E,^   A  +  C+I^l   v/2(l  +  6)    ^/   l+r  + 


5+2e 


^   yi-52  72+^(1  +  0  ]  72(1  +  6)    ^2(l  +  6)-(l  +  Oj 


(25) 


If  yxr  >>  I5  the  numerator  of  the  right  hand  side  of  equation 


(25)  Is  approximately  given  by 


/F 


greater,  using  crude  estimates,  than 
+1 


1 
/2l 


dC 


The  denominator  is 


1/2 


.1  ^^^   ^2+^(1-^^) 


1 
/27 


d? 


(26) 


1-5    ^3+C 


h.r        2  3/2/: 


so  that  the  right  hand  side  is  certainly  smaller  than  2 


3/2 


~T  >>    1  Is  incompatible  with  ^  >>  1.   If  j^  is  of  order 
unity,  but  bounded  away  from  zero,  we  see  that  the  numerator  of 
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the  right  hand  side  of  equation  (25)  is  of  order    —   while 
the  denominator  is  greater  than 


+1 


-1 


d5 


A-r2  TiT 


lf^(l+^) 


J2C 


Jl+l 


(27) 


In  this  event,  we  again  find  that  ^   can  not  be  much  larger  than 
unity.   We  conclude,  therefore,  that  yxT-  <<  1  is  a  necessary 

P  >>  1.   We  now  proceed  to  Investigate  this  limit. 


condition  for 


For  T  -*■  0  +  it  has  been  shown   that 
+  1 


d? 


-1 


-£^   ^l+£+T  ^ 


1  £n  32 

T 


1-5 


(28) 


Using  equation  (28)  and  t— r  <<  1  we  find 
+  1 


0         -,^1  +  ^     ^     n     1/2 

Zn    16 +    2e 

e 


2e 


1-C    ■^l+?+2c    Vl+C+^ 


:^  A+j;+2c  / 


d? 


/2(l+6) 


'l-C^  72(1  + 


(l+6)-(l+5) 


(29) 


There  are  four  regions  in  the  e-6  olane  in  which  t^  <<  1 

1  +  0 


These  are:   (1)  e  <<  1  and  6  <<  1;   (ii)  e  <<  1  and  6  =  0(1); 
(ill)  e  <<  1  and  6  >>  1;  (Iv)  e  >_  e   >  0,  6  >>  1  and  e  <<  6. 

For  case  (i),  we  easily  find,  making  use  of  equations  (A-1) 
and  (A-7)  of  the  Appendix, 


M-.l.i.n^  (30) 

^0 


^.l.l«„i6  ,3,) 


^=1,„16  (32) 

tj;     TT      6 


With  R  held  fixed,  we  see  that  this  limit  corresponds  to  i'^'^^- 

^2    H 
Moreover,  in  this  limit  -t—  '^  F-   I"  section  IV  we  shall  show 
A.     n 

that  this  limit  corresponds  to  the  case  where  the  column 
occupies  nearly  all  of  the  interior  of  the  container. 

Case  (ii)  also  corresponds  to  -^   large.   However,  -7—  =  0(1) 

1  ^0 

so  that  ijj^  =  0(1).   This  implies  that  A^  =  0(1)  while  A^  >>  1. 

The  compression  ratio   of  the  column  is,  therefore,  asymptoti- 
cally given  by  -r— 
^1 


^.1,— ^^^n±n <2       (34) 

1       f      dg         1 


^  A-?^   /2(l+6)-(l  +  5) 


/2(l+6)  ^    -,       p^„    r    ■>      1 

since  "ZIZIIZZIZZIIZ  -'■      ^°^   5    >   -1 

/2(l+6)-(l+C) 


12 


In  case  (ill)  we  obtain,  using  equation  (A-9)  of  the 
Appendix 


H   ,  1  ,„Lg)2  (l±6ii|  ^35^ 


^-2  (36) 

^.1,,  i  (37) 


T^  ^    1  •  (38) 


Moreover,  the  compression  ratio,  again  asymptotically  given  by 

•7—  ,  Is  asymptotically 
^1 

1^-2-1^      .  (39) 

A  A 

For  case  (Iv),  we  find  -r^  =  0(1).   For  e  >>  1,  -^  'v-  1, 


2-^  ^  el/2  ^^Q) 

^1 

the  column  contracts  to  a  nearly  circular  line  source.   In 

particular  ,„        ^ 

2£l/2  ^^   16  1+6 

ij;^       £n  I6e  ^^^^ 

1/2 

(12) 

(13) 
.n  166        •  "     ^''^ 
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R 

'X- 

27T    e"/^ 

K 

!?.n    166 

^2 

% 

2 

^0 

£n   I6e 

^ 

2 

Formulas  (30)  and  (31)  may  be  combined  to  give  an  asymp- 
totic estimate  of  &    in  terms  of  5-  and  e.   For  e  <  6  <  .001,  the 

n 

agreement  of  the  estimate  with  the  exact  value  of  6  is  better 
than  1% .      Similar  results  are  found  using  equations  (35)  and 
(36)  with  e  <  .01  and  6  >  1000,  and  using  equations  (4l)  and 
(42)  with  e  >  100  and  6  >  lOOe.   Referring  to  Fig. 7  ,  we  see 
that  the  above  regions  in  which  the  asymptotic  formulas  are 
accurate  cover  a  substantial  portion  of  the  figure.   The  level 


R 


The  asymptotic  expressions  for  -, —  and  -r—   determine  rela- 


tions  between  1+6  and  e   which  depend  only  on  the  ratio  5-.   Thus, 
in  case   (i)  and  in  case  (iv)  with  e  >>  1,  we  obtain 

1  +  6  ^  1  +  e^/^l6e^)^l-  «/")  (45) 

1  +  6  -.f^e^^  "/^)  (46) 

respectively.   These  formulas,  together  with  equations  (31)  and 
(42)  may  be  used  to  solve  for  e    and  1+6  as  functions  of  ii 
yielding  equations  (47)  and  (49).   Using  these  functions  of  i>^ 
we  obtain  equations  (48)  and  (50)  below.   Thus,  for  case  (i) 

1  +  6  '^.  1  +  I6e    °     ;     e  '^  l6e      °       (47) 


A,  -  ^,(^  -1) 

°  (48) 
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while  for  case  (Iv) 


1+6^1+  (l6e^)^l-^/"){i  f-  £n  (^  ?-)}^2R/H) 

IT  111  TT   ll;  ' 


A,  '^    ^^ 


A„  ^   i>^ 


(50) 


IT  ^) 


In  both  cases,  the  formulae  are  valid  for  the  limit  ^   ^0. 


single  parameter  tjj 
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IV.    Adlabatic  Compression 

We  now  consider  an  adiabatic  compression  of  the  plasma  by 
a  toroidal  discharge  at  the  outer  conductor.   We  shall  assume 
that  the  poloidal  field  (i.e.  the  field  in  the  plane  of  Fig.  2) 
vanishes  in  the  region  outside  the  outer  boundary.   With  this 
assumption,  if  I  is  the  magnitude  of  the  total  toroidal  current 
flowing  at  the  outer  boundary  and  B.  is  the  poloidal  field 
strength  at  the  interface,  then 

I  =  ^L  (51, 

where  L  is  the  length  of  the  perimeter  of  the  curve  bounding 
the  interface  in  the  plane  of  Fig.  2.   For  simplicity,  we  shall 
assume  that  compression  is  begun  with  a  column  filling  the 
container.   The  initial  strength  of  the  toroidal  field  trapped 
in  the  perfectly  conducting  column  will  be  given  by  B   .,  .   Its 
strength  at  any  other  time  will  be  indicated  by  B   , .   The  flux 
of  this  field  is  conserved  during  compression.   Therefore,  if  A 
is  the  cross-sectional  area  at  some  later  time 

B„^,  A  =  B*  ,  (4HR)      .  (52) 

col      col 

The  compression  ratio  of  the  column,  C,  is  then  given  by 

C  =  ^  .  (53) 

During  compression,  the  plasma  mass  is  also  conserved,  so  that 
the  density,  p,  initially  given  by  p  ,  satisfies 

pA  =  p*(4HR)  .  (54) 
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Assuming  the  compression  to  be  Isentroplc,  the  plasma  pressure, 
p,  and  temperature,  T,  are  computed  from 

P  ='^  P^  (55) 

T=a\^-^  (56) 

where  6<^  ,  c/   and  y   are  constants  {6(      =   where  p  =  pkT). 

It  Is  convenient.  In  place  of  the  pressure,  p,  to  define  a 
column  pressure,  p   .,  ,  for  the  plasma 

Pcol  =  P  "  27;  ^col    •  (57) 

The  initial  value  of  the  total  toroidal  current  flowing  at  the 

« 
outer  boundary,  I  ,  Is  given  by 


y        *^o^  col      \i  o^    2y     col' 


1 

'■0 


From  equation  (1) 


I  y^^col     L       .  (58) 


Since  p   .,  depends  only  on  the  initial  data  and  compression 
ratio  C   of  the  column  and  since  C  and  the  length,  L, 


ship  between  ^      and  I  must  be  found  numerically.   However,  for 
some  limiting  configurations  of  special  interest,  simple 
asymptotic  relations  are  readily  obtained. 
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If  the  height  of  the  column,  2A2,  satisfies  2A2/2H  %  1, 
the  pololdal  curvature  is  negligible  near  the  toroidal  mid- 
plane,  L/4H  ^  1,  C  ^  R/A-^  ,  and  ip^   ^   R-A.^.   For  y  =  2 ,  we  find 


~^  A~ 

I  ^1 

Using  equation  (^8)  we  obtain,  for  ^   /R  <<  1. 

R_  I 

I     '^       fi     .        7 


7-  ^  R 7  '     —  ^  I 7 

I      -| 1        ^O     — sr  -1 


(59) 


(60) 


From  equations  (3^),  (36),  (38)  and  (59),  we  find  that  maximum 
compression  ratios  for  a  column  with  high  eccentricity  are 
obtained  when  I/I   '\.  2. 

In  order  to  obtain  the  relationship  between  I  and  tj^   as 
the  column  cross-section  becomes  small,  we  first  obtain  the 
approximate  pololdal  field  at  the  toroidal  mldplane.   For  a 
rectangular  outer  boundary,  this  can  be  done  using  the  method 
of  Images.   In  particular,  if  H/R  >>  1 ,  we  may  neglect  end 
effects.   Replacing  the  column  by  a  line  source,  image  sources 
of  equal  strength  and  alternating  direction  are  added  at  y  =  J:2nR 
(see  Fig. 8  )  yielding  a  field  which  is  tangent  to  the  outer 
boundary  at  y=J:R.   On  the  y— axis  (the  toroidal  mldplane),  the 
pololdal  field  due  to  the  image  sources,  assuming  unit  current 
in  the  column  is 

^1  I    (_i)n     -2.y    1  (61) 

2^  [n=l      4n2R2_y2  J 
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The  associated  pololdal  flux  is  given  by 

£n  TT  (  -—2-  )       =  27  ^"  7  ^62) 

n=l     ^n 

In  obtaining  Eq .  (62)  the  true  dimensions  of  the  column  have 

been  neglected. 

If  the  column  has  a  cross— section  which  is  approximately 
circular  with  radius  A, ,  then  the  pololdal  flux  due  to  the 


column  is  p—  £n  v—  .   If  the  toroidal  current  flowing  in  the 
outer  conductor,  and  hence  the  column,  is  I,  (initially  given  by 


*    4H  /  ^ 

I  ^   —  -^2^  p   -,  )  ,  then  the  pololdal  flux  through  the  y— axis 

is  asymptotically  given  by 

2 
For  Y  =  2,  using  equation  (50),  and  A  "^   ttA,  ,  we  obtain 


iK) 

=v          ^      1               I 

I          „      o      R          n          R 

K    2  ^^  (I^) 

I                        ^0                *0 

(64) 
5-)       I"     %    % 
I 

Using  equations  (60)  and  (64)  we  obtain  the  following 
table  for  the  asymptotic  variation  of  the  indicated  thermo- 
dynamic parameters  of  the  system  in  terms  of  the  toroidal 
current  flowing  at  the  outer  boundary.   Figure  9  shows  the 
numerically  obtained  dependence  of  the  compression  ratio  on  this 
toroidal  current  when  5-  =  5 •9-   The  asymptotic  estimates 
of  the  compression  ratio  given  in  Table  II  are  also  shown. 
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Y    =    2 

p            T 

2^ 
P 

^2 
^1 

I 

rl     ^  (1+2R/H) 

I 

rl    -j2(l+2R/H) 
I 

(1       I    Ix    H    I 

I 

H/R    (I^^2 

i 

1 

TABLE  II 

As  seen  In  figure  9  the  compression  may  be  divided  into 

* 

three  stages.   Near  the  onset  of  compression,  I/I   less  than 

1.6,  the  column  maintains  a  shape  comparable  to  that  of  the 
outer  boundary.   For  the  Intermediate  compression  region,  with 
compression  ratio  between  2.5  and  10,  the  width  of  the  column 
remains  nearly  constant  but  its  height  diminishes  rapidly  as 
I/I   is  increased.   In  the  third  stage  of  compression  the  plasma 
cross— section  is  circular  and  very  high  compression  ratios  are 
obtained.   Plasma  cross— sections  representative  of  these  stages 
are  shown  in  figure  10 . 
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Appendix 

In  this  section,  we  derive  the  asymptotic  dependence  of 
the  numerator  of  equation  (25)-   The  relation 


+  1 


d^     

1-5^    A  +  5  +  T 


^—     -^^in^ 


(A-1) 


valid   for    t   ^-   0   +   vflll   be    used.      Employing   equation    (A-1),    the 


asymptotically    given   by 

I,    -v  ^  -^  iln{l6ii^  } 

^         /2(l+6)       /2  ^ 


(A-2) 


In  order  to  obtain  the  asymptotic  expansion  for  the  second 
Integral,  Ip,  we  Introduce  the  notation  i  =  2e,  e  =  y+t   • 


1/2 


y2(l  +  6)        J_^  A-K^      A  +  K  +  e      A  +  ?  +  e 
The  integral  in  equation  (A-3)  may  be  written  as 


(A-3) 


AL. 


y^/^/2-y   /y+i  /y+e 


0   ^s 


dy. 


J72 


y    /2-y   /y+e   /y+e 


(A-4) 
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where  s  Is  to  be  taken  "small"  as  prescribed  below.   If  s  is 
chosen  so  that  e/s^O,  we  find 


1/2 


M. 


1/2 


1/2 
y    /2-y   /y+e   /y+e 


^ 


1/2 


^2-y 
1/2 


1/2 
-]  /2 


(A-5) 


which  goes  to  zero  with 
go  to  zero  with  e. 


Therefore,  up  to  terms  which 


1/2 


I-.  "^ 


2   V2(l+6) 


dy 


y^^^/2^  Z^+f  /y+7 


1/2 


172 


dy. 


/2(l+6)   v^2~   „   y ^  /y+1  /y+e 


(A-6) 


The  change  of  variables  y  =  ev  yields 
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1//2 


s/e 


■2    /2(l+6)  Jq    v^/^(v+l)^/^(v+e/i)^/^ 


1//2       _^ 

/2TT+6T  ^Q   v^/2(v+l)^/^(v+e/e)^/2 


£/£ 


dv 


1//2 

/2TTT6T  ^_i      |v|^/2|^^;^,l/2|^^^/^|l/2 


dp 


l//2(l+6)   _^       - 

/  11-  e/e|   /2  J.i  yi-p^  ./l+p+  ^- 


/^ 


(A-?: 


z/Z 


If  ijj   ?^  0,  equation  (21)  implies  6  ?f  0,  so  that  e/e  7^  1.   In 
the  limit  1  +  6  ^«>,  equation  (A-1)  yields 


/2(l+6)   /2 


Jin  16(1  +  6) 


(A-8) 


In  this  limit,  the  denominator  of  equation  (25)  is  easily  seen 

to  be  asymptotic  to  —zzz   •      Thus, 
/l+T 


H^  1   ,n  (16)2  (llAV 


(A-9) 
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